Dark matter annihilation or unresolved astrophysical sources? 
Anisotropy probe of the origin of cosmic gamma-ray background 
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The origin of the cosmic gamma-ray background (CGB) is a longstanding mystery in high-energy 
astrophysics. Possible candidates include ordinary astrophysical objects such as unresolved blazars, 
as well as more exotic processes such as dark matter annihilation. While it would be difficult to 
distinguish them from the mean intensity data alone, one can use anisotropy data instead. We in- 
vestigate the CGB anisotropy both from unresolved blazars and dark matter annihilation (including 
contributions from dark matter substructures), and find that the angular power spectra from these 
sources are very different. We then focus on detectability of dark matter annihilation signals using 
the anisotropy data, by treating the unresolved blazar component as a known background. We find 
that the dark matter signature should be detectable in the angular power spectrum of the CGB 
from two-year all-sky observations with the Gamma Ray Large Area Space Telescope (GLAST), as 
long as the dark matter annihilation contributes to a reasonable fraction, e.g., > 0.3, of the CGB at 
around 10 GeV. We conclude that the anisotropy measurement of the CGB with GLAST should be 
a powerful tool for revealing the CGB origin, and potentially for the first detection of dark matter 
annihilation. 

PACS numbers: 95.35.+d; 95.85.Pw; 98.70.Vc 



I. INTRODUCTION 



What is the energy budget of the GeV sky? The origin 
of the cosmic gamma-ray background (CGB) in the GeV 
region, which was discovered by the Energetic Gamma 
Ray Experiment Telescope (EGRET) 0, Q, is a long- 
standing mystery. Unresolved blazars, a beamed pop- 
ulation of active galactic nuclei (AGNsLhave been the 
most popular explanation for the CGB d, 0, [E [E 0, B 
& Ell El however, even with the latest determi- 
nation of the gamma-ray luminosity function (GLF), it 
has been shown that only 25-50% of the CGB can be 
explained by unresolved blazars alone [lo| . Other astro- 
physical sources include clusters of galaxies, from which 
gamma rays are emitted by either hadronic collisions 
between shock-accelerated protons and the surrounding 
medium [l3l [Til [l5| or the inverse-Compton scattering 
of relativistic electrons off the cosmic microwave back- 
ground (cmb) photons [n, m m m, m nn m m. 

(Note, however, that such gamma rays have not been de- 
tected towards known clusters of galaxies yet [24| though 
there is some evidence [25|,[26j].) The CGB may be solely 
coming from these two unresolved astrophysical sources 
at the comparable level for each, or other sources may be 
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required to explain it. 

What could be the other sources of the CGB? R is 
energetically easy to produce gamma-ray photons in the 
right energy region by decay or annihilation of heavy (i.e., 
10 GeV or heavier) particles, the particles that have not 
been discovered yet. This possibility has been regarded 
as attractive because (i) more than 80% of matter in the 
universe is made of nonbaryonic dark matter, (ii) theories 
of supersymmetry predict that the lightest supersymmet- 
ric particles (e.g., neutralinos) are stable and are attrac- 
tive candidates for dark matter (in the mass range of 
GeV to TeV) , and (iii) these candidate dark matter par- 
ticles can annihilate into GeV gamma rays [13, [H, [2!| . 
The CGB flux from dark matter annihilation in cos- 
mological dark matter halos has been calculated in the 
GeV [3(| [U, HH as well as in the MeV energies jMSi. 
Successive works have shown that this mechanism can 
account for a large fraction of the CGB with a reason- 
able choice of dark matter parameters, with a significant 
boost of the signal from substructures |35| or minispikes 
around intermediate-mass black holes [36j (for the latter, 
see also Refs. [37], HU). In addition, with these boosts, 
we can also avoid an upper limit from the gamma-ray 
flux towards the Galactic Center [39[ . The excess of the 
CGB seen at around 3 GeV might be a signature of dark 
matter annihilation [401 ] . 

Ando and Komatsu have recently shown that the angu- 
lar power spectrum of anisotropy in the CGB may pro- 
vide a smoking gun discovery of the annihilating dark 
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matter [4l[ (hereafter AK06). In that study the authors 
focused on the annihilation signal from dark matter ha- 
los with smooth density profiles. Only semi-quantitative 
arguments were given for the other astrophysical sources 
and the effects of dark matter substructures. We there- 
fore investigate these two effects in detail in this paper, as 
they cannot be ignored if one wants to discuss whether 
anisotropy can really help detect the first signature of 
dark matter annihilation. 

This paper is organized as follows. In Sec. [Til we ex- 
plain what the CGB intensity averaged over all the di- 
rections looks like, for both the dark matter annihilation 
(Sec. Ill Ap and blazars (Sec. Ill B[) . We then turn our 
attention to the CGB anisotropy from dark matter sub- 
structure and blazars in Sees. [TTT] and IIVI respectively, 
where we present formulation and results of the angular 
power spectrum. Section[V]is the main part of this paper, 
devoted to discussion concerning anticipated anisotropy 
analysis in the presence of components from both dark 
matter annihilation and blazars. We study the case of 
other astrophysical sources and discuss the robustness of 
our results in Sec. IVII and also give conclusions in the 
same section. 



II. COSMIC GAMMA-RAY BACKGROUND: 
MEAN INTENSITY 



In this section we calculate the mean intensity (i.e., 
intensity averaged over the directions) of the CGB from 
both dark matter annihilation (with substructures taken 
into account) and blazars, and compare the characteris- 
tics of these two components. 



A. Dark matter annihilation 

We include the effect of dark matter substructures as 
follows: we assume that substructures consist of a num- 
ber of subhalos within a bigger host halo. These subhalos 
follow a certain mass function which is still unknown, but 
for our purpose we are only interested in quantities that 
are averaged over the mass function. If this mass func- 
tion is independent of the halo position as we assume, 
these averaged subhalos having the same gamma-ray lu- 
minosity would follow a smooth density profile of a host 
halo such as the one proposed by Navarro, Frenk, and 
White (NFW) [H, El with a halo concentration param- 
eter given in Ref. [4J|. The gamma-ray profile of a halo 
thus traces the dark matter density, rather than the den- 
sity squared which would be expected if dark matter dis- 
tribution were smooth [4l| . 

We define the number intensity, J/v, as the number 
of photons emitted per unit area, time, solid angle, and 
energy range. In a general cosmological context, it is 



given by 



EI N (h,E) = — dz 
4tt J 



P 7 ([l + z]E,z,hr) _ r(n 
H(z){l + z)* 6 



([l+z]E,z) 



where P 7 is the volume emissivity (energy of photons 
per unit volume, time, and energy range), H{z) 2 — 
Hq [(1 + z) 3 £l m + Qa] is the Hubble function in a flat uni- 
verse, and we assume the standard values for cosmologi- 
cal parameters, Ho = 100 h km s _1 Mpc -1 with h = 0.7, 
fl m = 0.3, and fl\ = 0.7. We specify a certain direction 
by a unit vector, n, position by a comoving distance vec- 
tor, r — rn, and time by a redshift, z (comoving distance, 
r, is also used interchangeably). The exponential factor 
reflects the effect of gamma-ray absorption due to pair 
production with the extragalactic background light; such 
an effect is negligible in the energy range of interest here. 

To evaluate the mean intensity, (In(E)), we need 
(Pj(E, z)). Let us define the gamma-ray spectrum per 
subhalo averaged over its mass function by M s \x{E), and 
the number of these subhalos within a parent halo of 
mass M by (N\M). Then, we obtain the mean volume 
emissivity as 



(Rr(E,z)) = (1 + z) 3 n sh (z)EM sh (E), 



(2) 



where n s h(^) is the mean comoving number density of 
subhalos given by 



M„ 



dM ^L(M,z)(N\M), (3) 



and dn/dM is the halo mass function for which we use 
the expression given in Ref. [45J| . The function, A4h (E), 
includes all the particle physics parameters such as the 
annihilation cross section, av, the dark matter mass m x , 
and the gamma-ray spectrum per annihilation dN^/dE. 
For the latter, we use a simple parameterization, i.e., 
dN ? /dE = (0.73/m, x )e~ 7 - 76E / m x /[{E/m^ + 0.00014], 
which is a good approximation for supersymmetric neu- 
tralino dark matter particles (30j . We parameterize the 
number of subhalos in each parent halo, which is also 
known as the Halo Occupation Distribution, as 



(4) 



If we ignore tidal destruction of subhalos entirely, we ob- 
tain a = 1, i.e., the number of subhalos is simply propor- 
tional to the mass of the parent halo. The tidal destruc- 
tion should also change the gamma-ray emission profiles 
in the parent halo, as it works more strongly at inner halo 
regions. However, we adopt the NFW profile for all of 
our calculations, because the profile change should exert 
only secondary effect to our conclusions as discussed in 

sec. mm 

Using Eqs. (fTJ) and 0, we obtain the mean intensity 
as follows: 



(In{E)) = / dr W([l + z]E,z), 



(5) 
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FIG. 1: The CGB spectrum from dark matter annihilation 
(dashed) and blazars with the best-fit LDDE GLF (dotted). 
Total intensity is shown by the solid curve, and the data points 
are from the EGRET data 0. 



where 

W(E, z) = ^-fi sh (z)Kh(E, z)e-« E ''\ (6) 

4-7T 

Figure [TJ shows the CGB spectrum from dark matter 
annihilation, where the particle mass is assumed to be 
m x = 100 GeV. We do not give a specific value of av or 
Mq. These parameters are degenerate, but they do not 
affect predictions of the angular power spectrum, as we 
see below. All we require here is that the predicted inten- 
sity becomes comparable to the observed CGB, and this 
can be done by adjusting these two parameters. Previous 
work which included dark matter substructures [35| has 
shown that this is indeed possible with a standard value 
of the annihilation cross section, av = 3 x 10 -26 cm 3 s" 1 , 
which gives the right amount of the dark matter density 
in the universe if dark matter was thermally produced in 
the early universe [27], . On the other hand, anisotropy 
depends only on m x and a. We shall therefore vary a 
and see how the results depend on a, while we fix the 
mass at 100 GeV throughout the paper. 



B. Blazars 

If a non-negligible fraction of the CGB flux comes 
from astrophysical sources such as blazars and clusters 
of galaxies, they inevitably give a background (noise) for 
the dark matter detection in the anisotropy signature. It 
is thus very important to evaluate the contribution from 



the unresolved point sources. We concentrate on blazars 
as an example. 

To calculate the mean CGB intensity from blazars one 
needs the GLF of blazars. We use the latest luminosity 
dependent density evolution (LDDE) model, which re- 
produces the observed GLF of the EGRET blazars bet- 
ter than a traditionally used pure luminosity evolution 
model flo| . As the LDDE GLF was originally given for 
the luminosity at 100 MeV, we need to generalize it to 
the other energies. We do this by specifying the spec- 
tral shape; here we assume it to be a power law with 
a spectral index of a 7 = 2.2 [l]. Then, the luminosity 
per unit energy range, C, is connected to the luminos- 
ity, L 7 (100 MeV) (= EC at 100 MeV) adopted in the 
previous GLF via the following simple relation: 



E C ni \ 1 

100 MeV J 



100 MeV ' 



( 7 ) 



The GLF is accordingly replaced with the one defined 
as the comoving number density per unit range in C, 
$>e(C, z), which is related to the original one through 



dC &e(C, z) — dL~ f p 1 (L 1 , z), 



(8) 



where we show the energy dependence of the new GLF 
explicitly by attaching subscript E. Note that p 7 on the 
right hand side is given by Eqs. (8) and (10) of Ref. [4(| . 
Using Eqs. (0 and (©, we can rewrite the luminosity and 
the GLF at any energies as long as the spectrum is kept 
to be a power law with the same index. 

The photon flux from the source with luminosity C at 
redshift z at energy E is given by 



(1 + z)C[{1 + z)E,z] 
4ndl(z) 



(9) 



where dL(z) is the luminosity distance out to a source at 
z. The flux sensitivity for point sources of the EGRET is 
F 7)Uro ~ 10" 7 cm" 2 s" 1 above 100 MeV and all the 
unresolved sources that give a flux below this threshold 
contribute to the CGB. The conversion from the differ- 
ential flux per energy, Te , to the integrated flux, F 1 , can 
easily be performed by integrating over energy above 100 
MeV and assuming the spectrum to be a power law with 
an index en-,. One obtains 



Te = (a-y - 1) 



E 



100 MeV 



(10) 



We use this equation and Eq. ([9]) to calculate the limiting 
source luminosity, C(J-E,iim, z), from F 7 ji m . 

We calculate the mean CGB intensity coming from un- 
resolved blazars whose gamma-ray flux is below J-E,\im 
from 



E(I N (E)} 



dz 



d 2 V 



i , dzdfl 



£(^ r E,lim,z) 



dC $ E {£,z) 

(11) 
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where we use z max = 5, and d?V/dzdVL is the comoving 
volume per unit redshift and unit solid angle ranges. We 
show in Fig. [T] the CGB spectrum calculated with the 
best-fitting LDDE GLF together with the EGRET data. 
The predictions fall below the EGRET data, accounting 
for only 25-50% of the observed CGB 0, Sf3|. 

This is presumably either because there is another class 
of objects which can contribute to the CGB by equally 
significant amount, or because the best-fitting LDDE 
model is underestimating the true GLF. 1 For the for- 
mer case, the particle acceleration in other astrophysical 
objects may also give a power-law CGB spectrum sim- 
ilar to that of blazars, and depending on its luminosity 
function, an unaccounted fraction of the CGB flux might 
be attributed to this population. Since the EGRET data 
may be explained by a power law component, such ad- 
ditional sources, plus blazars, might give almost full ac- 
count of the CGB in the GeV region. One such candi- 
date is galaxy clusters, in which protons or electrons are 
accelerated to relativistic energies by shock waves, and 
the GeV gamma rays are emitted by either pion produc- 
tion due to the proton-proton interactions [H, HI, HH 
or inverse-Compton scatte ring of relativistic electrons off 
CMB photons (IE IS EE Iial23, E3, S3, IS! ■ The latter 
possibility is also possible, although the GLF parameters 
that can account for 100% of the EGRET data are in- 
consistent with the X-ray data at the 4.4-er level. (See 
Sec. 2.2 of Ref. [46J for details.) 

GLAST is expected to detect 1000-10000 blazars as 
point sources [l(J, EH . We can thus improve on accuracy 
of the GLF significantly after GLAST with much better 
statistics. 



It is related to the spatial power spectrum of subhalos, 
P sh (k,z), through 



(I N (E)) 2 Q = J ^{W([l + z]E,z)} 2 P sh 



(14) 

where W(E,z) is again given by Eq. (JSJ), and P s h(k, z) 
may be divided into 1-halo (lh) and 2-halo (2h) terms: 



Psh(k) 


= P£( 

rOO 




J M m 



j2h/ 



dM ^L(mi) : 

dM \ n sh J 



\u(k\M)f 



Pl h (k) = 



dM ^L&l 

Mmin dM n sh 
x P lin (fc). 



b(M)u(k\M) 



(15) 



(16) 



(17) 



Each term means that we correlate two different points in 
one identical halo (lh) or two distinct halos (2h). Corre- 
spondingly, we also have 1-halo and 2-halo terms for the 
angular power spectrum, i.e., C; = C} h + Cf h . 

The 1-halo term is determined basically by the density 
profile of parent halos, p(r\M), as we assume that the 
subhalos in a parent halo distribute by following p(r\M). 
Here, u(k\M) is the Fourier transform of p(r\M)/M. The 
2-halo term, on the other hand, is proportional to the lin- 
ear matter power spectrum, Pu n (k) [49(, and depends on 
the halo bias, b(M) [E(J. Detailed derivations of these 
results are given in Appendix IA 11 We calculate the an- 
gular power spectrum from Eqs. (fT^ - (fT7|) with the NFW 
density profile of dark matter halos. 



III. COSMIC GAMMA-RAY BACKGROUND 
ANISOTROPY I: DARK MATTER 
ANNIHILATION 

The angular power spectrum of CGB anisotropy cal- 
culated by AK06 does not take into account the effect of 
dark matter substructure. In this section we extend their 
work by including substructures explicitly. 



A. Formulation 



The angular power spectrum, Ci, is given by 

Q = (\a lm \ 2 ), (12) 
f I N (E,h)-(I N (E)) 
ai m = J dil A (I N (E)) Y lmi n )- ( 13 ) 



1 This may instead be because of underestimated contamination 
of the Galactic foreground emission |4S| . 



B. Results 

Left panels of Fig. [5] show the results for two different 
Halo Occupation Distribution, a = 1 (top) and a = 0.7 
(bottom) [see Eq. flU)]. We have used M min = 1O 6 M 
in Eqs. p^|) and (fTT)) . The dependence on M m ; n is very 
weak, as the CGB flux from annihilation is dominated by 
massive halos which host many subhalos. 

The 1-halo term dominates at larger I's, or smaller an- 
gular scales, as expected. We find that the 1-halo term 
strongly depends on a: the smaller a, the larger the con- 
tribution from less massive halos. As smaller halos are 
dimmer, one needs to increase the source number in order 
to explain the observed CGB flux. Therefore, this makes 
the CGB more isotropic, resulting in a reduction of the 
1-halo term for a = 0.7. On the other hand, dependence 
of the 2-halo term on a is much weaker as the 2-halo term 
is essentially given by the linear matter power spectrum, 
Piin(k), times the average halo bias. This weak depen- 
dence comes from the bias factor that is in the integrand 
of Eq. (|17p . If we increase the contribution of low-mass 
halos by reducing a, then the 2-halo term decreases as 
such low-mass halos are less biased. 
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Multipole 1 

FIG. 2: Left: Angular power spectrum, 1(1 + l)C;/27r, of CGB from dark matter annihilation with substructures, (a) a. = 1 and 
(b) 0.7. The dotted and dashed lines show the contribution from the 1-halo and 2-halo terms, respectively, and the solid line 
is the total. Right: The same as left panels but with smooth dark matter distribution (i.e., no substructure; AK06), with the 
minimum mass of (c) 10 6 M Q and (d) 10 _6 M©. 



We use these two cases, a = 1 and 0.7, to bracket 
the uncertainty in a. With a = 1, we assume that all 
the substructures survive the tidal disruption, which is 
probably a very optimistic approximation. For a = 0.7, 
on the other hand, we allow some amount of disruption 
effect, and the angular power spectrum is totally domi- 
nated by the 2-halo term, another extreme situation. For 
even smaller values of a, the 1-halo term disappears but 
the 2-halo term decreases only slightly; thus we do not 
consider these cases any further. 

We also consider the case in which the smooth host 
halo component dominates the annihilation signal. (I.e., 
substructures are unimportant.) The gamma-ray emis- 
sion profile is proportional to the density squared, and 
this case was in fact carefully studied in AK06. (This 
case, however, requires additional assumption such that 
the density profile of the Milky Way is much shallower 
than NFW, in order not to violate the constraint from 
the Galactic Center [39].) The right panels of Fig.[2]show 
the angular power spectrum of smooth-halo cases with 
different values of the minimum mass, M m ; n . Unlike the 



subhalo-dominated case, the CGB flux is very sensitive 
to the choice of M m j n . This is because the halo concen- 
tration is larger for smaller halos [441 ] . and the gamma- 
ray luminosity from each halo is very sensitive to the 
concentration parameter. We again adopt two extreme 
cases here, one dominated by 1-halo term and the other 
by 2-halo term. The shape of the 1-halo term on the right 
panels is steeper than that on the left panels (subhalo- 
dominated case), as the signal is more concentrated at 
the central region. We see below that this tendency con- 
verges to C\ — const, when the source is infinitely small 
as expected for point-like astrophysical sources such as 
blazars. As for the 2-halo term, the shape is almost the 
same as the subhalo-dominated case. The normalization, 
however, is smaller, because the contribution from less 
massive halos is enhanced when the dark matter distri- 
bution is smooth (owing to large concentration of small 
halos), and less massive halos are less biased. 

As a final remark in this section, we discuss the case 
with other density profiles. Although the NFW profile is 
most widely used in the literature, our knowledge of the 
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density profile is still far from converged. In particular 
towards the central region, density might keep rising or 
converge at some constant value. Therefore, one might 
question about the impact of the different profiles on the 
angular power spectrum. We have already seen the ten- 
dency such that the steeper the gamma-ray emission pro- 
file becomes, the harder the 1-halo term of angular power 
spectrum gets. Exactly the same argument applies here 
as well. If the profile has a flat core within some radius, 
the 1-halo term should be flattened above some / that 
corresponds to the core radius, but this modification at 
such a small scale will not be detected with GLAST (see 
discussions about detectability in Sec. IV C[) . The flat- 
tened profile could also be caused by the tidal disruption 
of subhalos. This means that it controls both the shape 
(via the profile) and normalization (via the number dis- 
tribution, or a) of the 1-halo term; but the latter would 
be much more prominent. On the other hand, the 2-halo 
term would stay almost the same even if we changed the 
density profile, thus providing the guaranteed power spec- 
trum that is independent of the density profile adopted. 



IV. COSMIC GAMMA-RAY BACKGROUND 
ANISOTROPY II: BLAZARS 

A. Formulation 

When sources are point-like just as blazars, the angular 
power spectrum of the CGB is given by 



d = Cf + Cf, 



(18) 



where the first, Poisson term Cf , corresponds to the 1- 
halo term, and the second, correlation term Cf , to the 
2-halo term. The Poisson term represents the shot noise 
that does not depend on the multipole Z's, while the cor- 
relation term is due to the intrinsic spatial correlation of 
sources. These two terms are related to the spatial power 
spectrum through 



C 



dz 



d 2 V 



E 2 (I N (E)) 2 J " dzdfl 
x J dC <P E (£,z)F E (£,z) 2 , 



(19) 



I 



dzdtt lm I, r(z) ' 



E 2 {I N (E)) 2 
x J dC $ E {£, z)b B {£, z)T E {C, z 



(20) 



where the lower and upper bounds of integration over 
z and C are the same as in Eq. (fTTj) . Detailed deriva- 
tions are given in Appendix IA 21 Here the power spec- 
trum of blazars is approximated as Pb(/c; C\, C%) ~ 
&s(£i)fo B(£-2)P\in{k) , and the blazar bias, &b, represents 
how strongly blazars cluster compared with dark matter. 
(See also Sec. 2.1 of Ref. HI.) 



While the blazar bias, bs(C, z), is currently unknown, 
it will probably be measured directly from GLAST blazar 
catalog H|. (This measurement is, however, limited to 
the bias of resolved blazars, which can be different from 
the bias of unresolved blazars which contribute to the 
CGB anisotropy.) At the moment, one may estimate 
&s(£, z) from several approaches including the angular 
and spatial correlation analysis of optical quasars [Hit |52| 
and X-ray selected AGNs [H, H3, UK]. These results, 
however, are not consistent with each other, potentially 
due to some observations being biased by a limited field 
of view covered, or because there is something wrong in 
our understanding of the unified picture of the AGNs. 
In any case, a very wide range of the blazar bias is still 
allowed, b E < 5; see Sec. 3.2 of Ref. HI for a more 
detailed discussion. 



B. Results 

Figure[3]shows the angular power spectrum of the CGB 
from blazars predicted for EGRET. The dotted lines 
show the Poisson term [Eq. p9[) ] , the dashed lines show 
the correlation part [Eq. (f2"0"|) ] evaluated with bs = 1, and 
the solid lines show the total. While the blazar bias could 
perhaps vary from ~ 1 to 5, the Poisson term dominates 
the angular power spectrum at all multipoles if bs = 1. 
The dominance of Poisson term is due to a relatively 
small number of bright blazars just below EGRET's sen- 
sitivity. The Poisson term will decrease as we remove 
more fainter objects. We also remark that the C; does 
not depend on the gamma-ray energy, since we here as- 
sume all the blazars have the power-law spectrum with 
the same spectral index of a 7 = 2.2, and this energy de- 
pendence exactly cancels when we divide by the mean in- 
tensity squared {In(E)) 2 to obtain the normalized power 
spectrum C; [see Eqs. Jl9j and (|20]) ]. 

For GLAST, we choose the point source flux limit of 
2 x 10~ 9 cm~ 2 s- 1 (E > 100 MeV), ~ 50 times bet- 
ter than EGRET, which is expected to be achieved after 
two years of all-sky survey observations of sources with 
a spectral index of 2 H|. Our predictions for Ci from 
GLAST data are shown in Fig. [4] As GLAST can detect 
and remove more fainter objects than EGRET, the Pois- 
son term is greatly reduced while the correlation part is 
almost unchanged. If the blazar bias is larger than 1, 
the correlation part would dominate the angular power 
spectrum at low Vs, which would allow us to measure the 
average bias of unresolved blazars. 

We also show the correlation part of the angular power 
spectrum using a bias model which was inferred from the 
optical quasar observations [U : 



b Q ( z ) = 0.53 + 0.289(1 + zf 



(21) 



If the unification picture of the AGNs is correct, then 
it may be natural to set bs = 1>q(z). The results from 
this calculation are shown as the dot-dashed curves in 



7 




1 10 10 2 10 3 

Multipole 1 



FIG. 3: Angular power spectrum of the CGB from unresolved 
blazars expected from the EGRET data. Contributions from 
Poisson term, Cf , and the correlation term, Cf with 6b = 1 
(pB = bcj(z)), are shown by the dotted and dashed (dot- 
dashed) curves, respectively. The total contribution is shown 
as the solid curve for bs = 1 



Figs. [3] and |4j We find that these results are quite sim- 
ilar to the case of b B = 1- This is because at low red- 
shift, z < 0.5, the quasar bias is close to 1, and the main 
contribution to the CGB from blazars comes also from 
relatively low-redshift range. Once again, we note that 
the quasar bias [Eq. (f2"Tj) ] is significantly different from 
the bias inferred from the X-ray AGN observation, which 
indicated stronger clustering [EH, [H, HH . Therefore, one 
should keep in mind that a wide range of the blazar bias, 
possibly up to ~ 5, is still allowed. Hereafter, we adopt 
b B = 1 as our canonical model, and note that Cf simply 
scales as b 2 R . 



V. DISTINGUISHING DARK MATTER 
ANNIHILATION AND BLAZARS 

The main goal in this paper is to study how to dis- 
tinguish CGB anisotropics from dark matter annihila- 
tion and from blazars. The current uncertainty in the 
blazar bias would be the source of systematic errors, but 
this can be reduced significantly by several approaches, 
such as the upgraded and converged bias estimations of 
AGNs from the other wavebands, direct measurement 
of the blazar bias from the detected point sources by 
GLAST Hi| , and the CGB anisotropy at different ener- 
gies where the contribution from dark matter annihila- 
tion is likely to be small. 



A. Formulation for the two-component case 
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FIG. 4: The same as Fig. [3] but for the CGB anisotropy ex- 
pected from GLAST data. 



The total CGB intensity is the sum of dark matter 
annihilation and blazars: 



I CGB (E,n) 
</cgb(£)> 



I B (E,n) 



I D (E,n), 
(Id(E)), 



(22) 
(23) 



where the subscripts B and D denote blazar and dark 
matter components, respectively. The expansion coeffi- 
cients of the spherical harmonics are given by 



CGB 
Hm 



dflf 



I GGB (E,n)~ (I GGB (E)) 

{Icgb(E)} 
6I B (E,n) + 8I D (E,n) 



YL(n) 



(Icgb(E)) 



(24) 



where SI b .d = Ib,d - (Ib,d), Ib,d = (Ib,d)/ (Icgb)- 
These f B and fo are the fraction of contribution from the 
blazars and dark matter annihilation to the total CGB 
flux, and we have the relation f B + /d = 1- There- 
fore, af r ^ D is defined as the coefficient of the spherical 
harmonic expansion if each component is the only con- 
stituent of the CGB flux, the same definition as in the 
previous sections or of AK06. The total angular power 
C, CGB = (|a^ B | 2 ), is therefore written as 



spectrum, 
£«CGB 



f B Cl. B + f B Ci,D 



^Ib/dCi^d, 



(25) 
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where Ci >E and Ci iE are the angular power spectrum of 
the CGB from blazars (Sec. IIV[) and dark matter anni- 
hilation fSec. IIIII and AK06), respectively, and Ci.bd = 
( a hn a im) ^ S a cross correlation term. This cross correla- 
tion term is derived in Appendix [Bl and is again divided 
into 1-halo and 2-halo terms, i.e., 



-<2h 



(~iVh i r-i2 

l,BD — ^LBD T ^l.BDi 



(26) 



where each term is given by 

W([l + z]E,z) 



CIbd = I dr 



x F E {C,z 



{I B {E)){I D {E)) 

(N\M(C)) (I 



J dC $ E (C,z) 



C 2h 



BD 



dr 



W([l + z]E,z) 



M[L] , (27) 



(I B (E)){I D (E)) 



x T E (C,z)b B (C,z) / dM 



x b(M,z)u ( l -,z;M) Pun i-,z 



for the subhalo-dominated annihilation, and 
W([l + z]E,z) 



dn(M, z) (N\M) 
dM n sh (z) 

(28) 



Cl h BD = I dr 



E(I B {E))(I D (E)) 



J dC $> E (C,z) 



c, 



2h 

l.BD 



xF E (£,z)^-v(-,z;M[£]), ,2!)) 



dr 



^mPc \r 

W([l + z]E,z) 



E(I B (E))(I D (E)) 



^jdC^Cz) 



x F E (C,z)b B {C,z) / dM 



dn(M,z) M 



dM n m p c 

(30) 



. hi V. [ 1 -,z-,m\ P lin ( l -,z 



for the host-halo-dominated annihilation; v(k, z) is the 
Fourier transform of p 2 (r\M)/Mil m p c . Here we note 
that the function W(E, z) in the latter expressions 
[Eqs. ([23) and 1(50])] is different from Eq. ©, but is given 
by Eq. (5) of AK06. In order to evaluate the 1-halo term, 
we need a relation between blazar luminosity C and its 
host mass M, for which we use Eq. (22) of Ref. [46f . 



B. Anisotropy due to dark matter annihilation in 
the two-component case 

Since our main thrust here is how to detect the dark 
matter annihilation signature out of the CGB in the pres- 
ence of more common (and plausibly known) blazar com- 
ponent, we focus our attention only on the energy of 10 
GeV, a typical energy of gamma rays expected from the 
annihilation of dark matter particles with a mass of ~ 100 
GeV. From now on we treat the blazar contribution as 
the "background noise." A signal and background of the 



angular power spectrum is, respectively, 
C? = f D Ci,D + 2f D {l-f D )C i 



LBD, 



LB, 



(31) 
(32) 



and we assume that the background is very well known. 
This is a reasonable assumption, provided that we can 
model the GLF and bias of unresolved sources from those 
of resolved (detected) sources in the GLAST data. One 
can also calibrate the background (blazar) component 
using the angular power spectrum of the CGB at lower 
energies such as 100 MeV, where contribution from the 
dark matter (neutralinos) annihilation may be ignored. 

We use several LDDE parameter sets for the blazar 
GLF, which are listed in Table |TJ The faint end of the 
GLF, which makes the largest contribution to the CGB, 
is given by p~ i oc kL^ 71 , where 71 is the faint-end slope 
and k is the overall normalization [1 01 ] . We vary these pa- 
rameters such that the blazars explain a certain fraction, 
/| GRET , of the CGB intensity measured by the EGRET 
at 10 GeV. We call these models LDDE10, LDDE30, 
LDDE50, and LDDE70 for /| GRET = 0.1, 0.3, 0.5, and 
0.7, respectively. We note that the best-fitting LDDE 
model, (71, k) = (1.19, 5.11 x 10~ 6 ), explains only - 20% 
of the CGB flux at 10 GeV, and therefore, one should 
keep in mind that this list includes models that are some- 
what disfavored from the redshift and luminosity distri- 
butions of the EGRET blazars. In the two-component 
treatment, the dark matter fraction is obviously obtained 
by /d = 1 - Ib- 

As for the GLAST data, the contribution to the CGB 
from blazars will be greatly reduced as GLAST can re- 
solve and remove more blazars. On the other hand, 
the contribution from dark matter annihilation will not 
change because it is unlikely that we can detect in- 
dividual dark matter halos via annihilation even with 
GLAST [H [H. Therefore, the total CGB intensity, 
(-^CGb), that would be measured by GLAST would be 
smaller than that measured by EGRET, while keep- 
ing intensity from dark matter annihilation unchanged. 
The expected fractional contributions to the CGB in the 
GLAST data from blazars or dark matter annihilation 
are also shown in Table HI 

In Figs. [5] and [6) we show the angular power spectrum 
of dark matter annihilation (i.e., Cf\ solid) as well as 
of background (Cf; dotted) that would be measured by 
GLAST, for the blazar models shown in Table U We 
assume that the dark matter signal is dominated by sub- 
structures. We also show the expected 1-er error bands 
of Ci . Note that the error bands do include contribu- 
tions from the blazar power. (See Sec. IV CI for the error 
estimation.) The subhalo distribution in each halo is as- 
sumed to be (N\M) oc M (Fig. [5]) and oc M - 7 (Fig. O. 
The shape of the angular power spectrum is quite dif- 
ferent between dark matter annihilation and blazars, a 
characteristic that should be useful for a smoking gun 
detection of dark matter annihilation. 

Even though the dark matter contribution to the CGB 
is relatively small for a = 0.7, the cross correlation term 
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TABLE I: Blazar LDDE GLF models adopted in the two- 
component (dark matter annihilation and blazars) study of 
the CGB anisotropy. Values of 71 (the faint-end slope of GLF) 
and re (the overall normalization of GLF) specify each model, 
and the output is given in terms of the fractional contributions 
to the CGB from blazars, fs, and dark matter annihilation, 
fn, for EGRET or GLAST. 



Model (7i,10 6 re) 


rEGRET 
JB 


(EGRET 

Id 


/•GLAST 
JB 


/•GLAST 
JD 


LDDE10 (1.05,6.33) 


0.1 


0.9 


0.03 


0.97 


LDDE30 (1.23,4.72) 


0.3 


0.7 


0.20 


0.80 


LDDE50 (1.29,4.11) 


0.5 


0.5 


0.39 


0.61 


LDDE70 (1.33,3.69) 


0.7 


0.3 


0.61 


0.39 



(due mainly to 2-halo term) is still reasonably large; we 
find that the 1-halo term of the cross correlation is always 
negligibly small as long as an empirical luminosity-host 
mass relation is used (Eq. (22) of Ref. [46|). Although 
the shape of the cross correlation is similar to the pure 
blazar correlation term (because both are proportional to 
the linear power spectrum Pr m ), its normalization gives 
us useful clue of the source because the normalization 
of the pure blazar term should be known to reasonable 
accuracy in advance. 

We then repeat the same analysis for the host-halo- 
dominated annihilation (no substructures). This is a 
straightforward extension of the study in AK06, to which 
the background is added. We show in Figs. [7] and [5] the 
power spectra for M m ; n = 10 6 Af Q and 10 _6 M Q , respec- 
tively. The smaller minimum mass is motivated by the 
free-streaming scale of neutralinos, 1O~ 6 M [57l. IBa. l59]| . 
while the larger one represents an extreme case in which 
halos smaller than a million solar masses (M < 10 6 M Q ) 
have been tidally disrupted. (Note, however, that 
there still remains larg e allowed range for the former 
case, 10 _12 -0.1Mq [59!]). Compared with the subhalo- 
dominated case, the anisotropy signature is typically 
smaller, but the general tendency is almost the same, 
justifying qualitative arguments regarding substructures 
given in AK06. 

C. Can GLAST detect dark matter annihilation? 

We use the standard procedure to calculate the pro- 
jected 1-er error (binned over A^) on the extracted power 
spectrum of the CGB from dark matter annihilation: 

f ^/S( c;+c '' + l)' (33) 

where / s k y = f2 s ky/47r is the fraction of the sky cov- 
ered by GLAST, Al is the bin width (which we shall 
take to be Al — 0.51), CV = ^ sky A^ total /A^Q GB is the 
power spectrum of photon noise, and A'cgb and A^otai 
are the photon numbers of the CGB and total (CGB plus 



other backgrounds), respectively, expected from the re- 
gion r2 s ky, and Wi is the window function of a Gaussian 
point spread function, Wi — exp(— Z 2 er 2 /2). Note that 
this formula assumes that CGB anisotropy obeys Gaus- 
sian statistics. We take the following specifications for 
GLAST: the field of view is f7f ov = 2.4 sr, the angular res- 
olution is a;, = 0.115°, and the effective area is A e g = 10 4 
cm 2 , both evaluated at E — 10 GeV [56j. We assume a 
two-year all-sky survey (T = 2 yr) , which corresponds to 
mean exposure of A c gt c s — A c gTflf ov /4:ir = 1.2 x 10 11 
cm 2 s, towards each point in the sky. 

Equation clearly shows that the "astrophysical 
background noise" from blazars, Cf, contributes to the 
error budget. On the other hand, the background that 
contributes to Cat includes detector noise and Galactic 
gamma-ray radiation, which we call foreground. The 
detector noise is very small for GLAST, about 5% of 
the CGB flux above 100 MeV; this decreases signifi- 
cantly as the gamma-ray energy and the detector noise 
is safely assumed to be negligible. The Galactic fore- 
ground is more difficult to estimate, but according to 
recent calculations using a numerical code of the Galac- 
tic cosmic-ray propagation, its intensity at 10 GeV is 
estimated to be E 2 lN tga \ ~ 10~ 7 GeV cm~ 2 s _1 sr -1 
at high Galactic latitudes, \b\ > 20°, which is an or- 
der of magnitude smaller than the observed CGB inten- 
sity (see Fig. 4 of Ref. [60]). Below this latitude, the 
Galactic foreground entirely dominates the gamma-ray 
flux; thus, we do not consider \b\ < 20° further. We 
also note that the Galactic foreground may contribute to 
the C\ due to spatial fluctuation of the sources, but we 
neglect this effect in the following arguments. We pos- 
tulate that this is a reasonable approximation because 
its flux at the large Galactic latitude is very small com- 
pared with the CGB, which suppresses the contribution 
to the total C/ as well; including this effect requires a de- 
tailed modeling of the Galactic cosmic-ray propagation, 
which is beyond the scope of the present study. The 
fraction of the sky relevant to our analysis is therefore 
/sky = 1 — cos 70° = 0.66, which is reasonably close to 1. 
We do not lose much by this Galactic cut. The number 
of CGB photons expected from this region and for T = 2 
yr is Ncgb =s EI CG BA cS t cS n sky = 10 5 (£/10 GeV)" 1 , 
while the corresponding photon count from foreground is 
~ 10 4 — negligible compared with that of the CGB. We 
finally obtain C N ~ n sky /N CGB = 8 x 10~ 5 (£;/10 GeV) 
sr. Note that this number is based upon the CGB inten- 
sity measured from EGRET, and it will be smaller for 
GLAST because GLAST can detect and remove more 
fainter blazars. While this reduction is eventually found 
to give no significant effect, we take more precise ap- 
proach in the actual error estimation. Here, just in order 
to provide a rough idea of the size of each quantity, we 
used the CGB flux obtained with EGRET. 

We show our predictions with the error bars calculated 
from Eq. ([55)1 in Figs.JSHHl We find that, for any models 
we have considered here, one should be able to detect 
the angular power spectrum from dark matter annihila- 
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FIG. 5: Angular power spectrum of the CGB from dark matter annihilation (/cG.d; dashed), blazars (f%Ci,B\ dotted), and 
cross correlation ^/b/dCbd; dot-dashed) that would be measured by GLAST at E — 10 GeV, for various models of the blazar 
GLF and various fractions of dark matter contribution to the CGB, fn (Table [TJ) . The adopted dark matter mass is 100 GeV 
and the gamma-ray emission is assumed to be dominated by the substructure. The total signal, Cf = fr>Ci,D + ^Jb/dC^bd, 
is shown as the solid curve, while the corresponding GLAST errors (<5(7f; for two years) are indicated as boxes. The signal is 
to be detected if it is larger than the size of errors (Cf > SCf). The subhalo distribution in a halo of mass M is assumed to be 
(N\M) oc M. 



tion with GLAST in two years of operation, as long as 
the dark matter contribution to the mean CGB flux is 
greater than 30% at some energy within the GLAST en- 
ergy window. This statement is independent of the den- 
sity profile adopted as it is derived with the guaranteed 
power spectrum, the 2-halo term. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we have presented detailed calculations 
of the angular power spectrum of CGB anisotropy from 
dark matter annihilation in cosmological dark matter ha- 
los as well as from unresolved blazars. The power spec- 
trum of dark matter annihilation from smooth NFW 
halos (he-, no substructures) has been calculated by 
AK06 [4l( , and the power spectrum of resolved (detected) 
blazars has been calculated by Ref. [46[ . Our work builds 



on and extends these results by taking into account the 
effects of dark matter substructures explicitly, by means 
of the Halo Occupation Distribution of subhalos. These 
calculations should provide a useful benchmark for the 
angular power spectrum of CGB that would be measured 
by GLAST. 

Our results are very encouraging: one should be able 
to detect the angular power spectrum from dark mat- 
ter annihilation with GLAST, whether dark matter halos 
are smooth or clumpy, as long as the dark matter con- 
tribution to the mean CGB flux is greater than 30% at 
some energy within the GLAST energy window. This is a 
rather modest requirement given the fact that unresolved 
blazars appear to contribute to the mean CGB only at the 
25-50% level As far as the mean CGB is concerned 
it has been pointed out that subhalos are necessary in 
order for dark matter annihilation to make a significant 
contribution without violating stringent constraints from 
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FIG. 6: The same as Fig. El but for (N\M) oc M° 



the Galactic Center 35, 36, 39]. Thus, our current "best" 
predictions are either Fig. [5] or Fig. O depending on the 
degree to which tidal disruption of subhalos is important. 

The "background noise" for dark matter annihilation 
searches includes the blazar anisotropy, detector noise, 
and the Galactic foreground. The blazar anisotropy 
should be well calibrated with the CGB anisotropy at 
lower energies (where dark matter signal is probably 
unimportant), analysis of AGNs selected with other 
wavebands, and/or the power spectrum of resolved point 
sources. The detector noise is always negligible. The 
Galactic foreground contamination is serious near the 
plane, \b\ < 20°, whereas its flux is found to be at most 
~ 10% of the CGB flux at high Galactic latitudes. The 
anisotropy due to the foreground is thus accordingly re- 
duced. We also note that dark matter annihilation in the 
Galactic halo may also give a contribution to the CGB 
flux comparable to that from the cosmological halos (e.g., 
Ref. [H,|6l|). Adding this component would increase the 
predicted anisotropy from dark matter annihilation. 

GLAST covers energy spectrum up to 300 GeV, which 
might enable us to probe a line signature in the CGB 
due to direct annihilation of dark matter particles into 
two photons (or one photon plus one Z° boson). In ad- 



dition to quite robust spectral feature, one could also use 
the CGB anisotropy as a consistency check. When the 
spectrum has a feature as the case of lines, the angu- 
lar power spectrum becomes larger [39l l62j , which makes 
this method even more promising. 

Finally, we comment on how our conclusions might 
be affected by the other astrophysical sources (of either 
known or unknown species) that we have not considered. 
If the gamma-ray emitter is a point source and follows 
the distribution of dark matter halos, then the number of 
such sources contributing to the CGB is the only param- 
eter that determines the amplitude of the power spec- 
trum. Whatever the sources are, they should be dimmer 
than blazars on average; we must have seen them in the 
EGRET data otherwise. If the number of such sources 
is much larger than blazars, then the Poisson term, Cf, 
is highly suppressed and is unimportant. If these sources 
are rare, forming only in large halos, the power spectrum 
(both 1-halo and 2-halo terms) can be large, but it is ex- 
tremely difficult for them to give a dominant contribution 
to the CGB. 

If these sources are spatially extended, e.g., clusters of 
galaxies, then the shape of the angular power spectrum 
depends on their gamma-ray profile. How the galaxy 
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FIG. 7: The same as Fig. [5] but in the case of the host-halo-dominated annihilation (no substructures; AK06). The minimum 



mass of dark matter halos is assumed to be 10 M, 



0- 



clusters are extended in gamma rays depends on the 
emission mechanism, which is uncertain. (No direct de- 
tection of gamma rays towards known clusters has been 
made so far.) Incidentally one can make the most con- 
servative (or pessimistic) prediction as to how much the 
unresolved clusters would contribute to CGB anisotropy, 
by assuming that the galaxy clusters are point sources. 
We have performed the same calculations for galaxy clus- 
ters for both the proton-proton collision [3| and inverse- 
Compton models [T3|. Our results again suggest that 
the dark matter component should be detectable signifi- 
cantly, as long as its contribution to the EGRET CGB is 
more than ~ 30%. Angular distribution of background 
radiation from galaxy clusters has been investigated also 
in previous papers, for radio waveband [HI, [63| as well as 
for gamma rays [TH, [l9[ . 

Based upon these results, we conclude that anisotropy 
in the CGB that would be measured by GLAST has to 
be analyzed in search for signatures of dark matter anni- 
hilation. If detected, it should provide us with the first 
(albeit indirect) evidence for emission from dark matter 
particles, which would shed light on the nature of the 
still-mysterious dark matter in the universe. 
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APPENDIX A: RELATION BETWEEN 
ANGULAR AND SPATIAL POWER SPECTRUM 

1. Dark matter substructure 

In this subsection, we derive Eq. (fl"4|) following the halo 
model approach [64 . 65] . Since the gamma-ray intensity 
due to annihilation dominated by subhalos depends on 
the density, which can be written as a superposition of 
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FIG. 8: The same as Fig. [7] but with minimum halo mass of 10 6 Mg 



densities in different halos, we obtain the following ex- 
pression for the volume emissivity: 

P 7 (E, x) = I dM 1 d 3 x 1 5 d(Mi - Mi)6 3 D (xi - a*) 

x (jV|M)u(as - xi\M 1 )EJ\f sh (E) (Al) 

where <5^ is the iV-dimensional delta function, and x 
represents comoving coordinate. As in Ref. [64[, the en- 
semble average of the sum over delta functions is simply 
the seed density: 

l^5Ux 1 -x t )8 D {M 1 -M l )\^ d -^, (A2) 



and using this into Eqs. (|A1[) and ([T|), we recovers Eq. ([5]) 
since J d?X\u(x — x\\M\) = 1. 

We then define the two-point correlation function of 
subhalos as 



£ sh (a: - y) 



(SR^SP^y)) 
(P 7 (|x|)>(P 7 (|y|)>- 



(A3) 



where <5P 7 (a:) = P-y(x) — 5P y (\x\). Now we evaluate 
(P~ ( (x)P 1 (y)). The ensemble average of the product of 
seed densities is written as follows 1641] : 



dn(Mi) dn(M 2 
dM dM 



l + ^(x 1 -x 2 -M ll M 2 ) 



dn(Mi) 
dM 



5%{x 2 -xi)fe(M 2 -Mi), 



(A4) 



r 



where %i N ^ is the TV-point correlation function of the seed. The first term represents the two- halo contribution, i.e., 
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two points considered, X\ and ai 2 , are in two distinct ha- 
los, while the second are the one-halo contribution where 
these two points are in the same halo. By substitut- 
ing this expression into Eq. (|A1[) , and using the relation 
(SP y {x)6P 7 (y)) = (P 7 (z)P 7 (y))-<P 7 (M))<P 7 (|y|)} and 
Eq. (©, we obtain 

dM \ n sh / 
x u(x — x\\M\)u(y — xi\M\) 



d 3 x 1 d 3 x 2 dM 1 dM 2 



dn{M x ) dn(M 2 ) 
dM dM 



(JV|Mi)(JV|M 2 ) , , , 

x — — u(x - a;i Mi) 

x u(y-x 2 \M 2 )^(x 1 -x 2 ;M 1 ,M 2 ). 

(A5) 

We define the power spectrum of subhalos P s h(k) as 
Fourier transforms of £ s h(r)- Remembering that convo- 
lution in the real space corresponds to a simple product 
in the Fourier space, we obtain for the expression for 



Psh = P 



Ih 
sh 



Pju 1 as follows: 



P^(k) = J dM 
Pl h (k) = 



dn(M) f(N\M) 



dM 



dM \ n sh 

dn(M) (N\M) 
dM fi s h 



,T\\2 



\u{k\M)\ 



(A6) 



b{M)\u{k\M)\ 



X Pin(fc), 



(A7) 



(2) 

where we used an approximation as Q (f", Mi, M 2 ) « 
b(Mi)&(M 2 )£iin(j") by introducing the bias factor b(M) 
and linear matter correlation function £i; n . These ex- 
pressions are identical to Eqs. (fH)|) and (p~7|) . 

Finally we derive the relation between angular power 
spectrum C; and spatial power spectrum P a h(k) in the 
following. We start with the angular correlation function 
that is defined by 



(I N {E)) 2 C(0) = (8I N (n 1 ,E)6I N (h 2 ,E)), 



(A8) 



where cos 9 = fi\-fi 2 . Using this definition together with 
Eqs. (J!]), ©, @, and ([AH]) , we get 

(/jv) 2 ^) = / A-xdr, W(zi)W(z 2 )£ sh (n -r 2 |z!,z 2 ), 



(A-9) 

where zi and z 2 are redshifts corresponding to T\ and r 2 , 
respectively, and we suppressed energy indices for sim- 
plicity. To further simplify, we use small separation ap- 
proximation following Ref. [6(| [we use dr\dr 2 = drdr], 
with definitions of r = (r*i + r 2 )/2,rj = r 2 — r{\, and 
arrive at 

(In) 2 C(8) = f drd v W(z) 2 U(vr + rO\z). (A10) 



The angular power spectrum is related to the correlation 
function through 



C,= d 



j2q -ue 



C(6), 



(All) 



for small scales, or for large multipoles. Then, the Fourier 
transformation of the relevant quantity simplifies as fol- 
lows: 



d 2 9 e- u e 
= I d~V i 



2a „-u-e 



d 3 k 



drj 



(2tt) 2 



P s h(fc, z)e 



d 2 / dr) 



*^P sh (fc,z)e^»V e -(^-') 



J dk {l d 2 k± P sh yk 2 + k 2 ± , z) 5 D (k {{ )5 2 D (rk x - I) 

(A12) 



\-P s h ( k = -,z 



In the second equality, we decomposed the wavenumber 
k by the components parallel and perpendicular to r, 
i.e., k = ft 1 1 + fej_, and used d 3 k = dk^d 2 k±. Therefore, 
together with Eqs. (jAlOl) and (jAl 1[) . we arrive at our 
demanded relation, Eq. (Q2 



2. Blazars 

In this subsection, we derive Eqs. (fTT))) and (|2^|) fol- 
lowing discussions in Ref. (see also, Ref. [67j). The 
Poisson noise Cf is obtained by setting 8 = for C{6) 
[Eq. flM}], and it is given by Eq. (58.14) of Ref. [gS| [note 
with Eq. (58.12) that the definition of the angular power 
spectrum is slightly different] , which is 

E 2 (I N (E)) 2 Cf = Jdz^JdC 3(4 z)T E (C,z) 2 , 

(A13) 



in our notation. It is equivalent to Eq. |T| 

The correlation term of the angular power spectrum 
Cf , on the other hand, is obtained with the Fourier 
transformation of the angular correlation function for 
6 ^ (Eq. (58.6) of Ref. Hi) 



E 2 {I N (E)) 2 C(6) 



1 



16tt 2 



dz 



d 2 V 



dzdn (1 + z) 2 r(z) 2 J_ 



du 



x £ B {ur + r(z)66,z) 



dC £$ B (£,z) 



(A14) 



where £s(r, z) is the two-point correlation function of 
blazars; here we suppressed the luminosity index £ but 
note that this dependence should appear in general. 
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Then, using the Fourier transformation that is similar to 
Eq. (|A12|) . we obtain the correlation part of the angular 
power as 



E*(I N (E)) 2 Cf 

d 2 9 e~ a - 



= / dz 



d 2 V 
dzdQ 



'E 2 (I N {E)) 2 C(6) 
1 



-1 2 



J dC C<S> E {C,z) 



4ir(l + z)r(z) 2 

2 



Pe 



r(z 



dz - — —Pb 



dzdQ " \t{z) 
J dC <S> E {C,z)T E {C,z 



(A15) 



where in the last equality we used the relation r{z) = (1 + 
zy 1 d L (z) andEq. ©. This is very similar to Eq. f20]). It 
is obvious that if we introduced the bias parameter from 
the beginning, as £ B (r;£i,£ 2 ) = b B {£i)b sO^dinO"), 
then we would have obtained exactly the same result as 
Eq. ®. 



APPENDIX B: CROSS CORRELATION 
BETWEEN DARK MATTER ANNIHILATION 
AND BLAZARS 

We here derive the formulation for the angular cross 
correlation between blazars and dark matter annihila- 
tion, Ci^bd, in the subhalo-dominated case, Eqs. ([27)) 
and l|28p. and in the host-halo-dominated case, Eqs. (j2"5)) 
and (|30p . respectively in the following subsections. Once 
again, we follow the halo model approach highlighted in 



Refs. |64|,l65j. 

1. Subhalo-dominated case 

In the halo model approach, the blazar intensity can 
be written as 

EI B {E,n) = — J d 3 x J dC J2 6 d( x - x i) 

x S D {£- Ci)^ E {C,z), (Bl) 

and the intensity from dark matter subhalos are Eq. ([T]) 
with P 7 given by Eq. (|A1[) . Therefore, the procedure of 
obtaining (IbId) is quite similar to that in Appendix IA II 
Repeating arguments there, we arrive at Eqs. (|27|) and 

2. Host-halo-dominated case 

In the case of host-halo-dominated annihilation, the 
intensity is proportional to density squared: 

EI D (E 7 h) = J dr ES 2 {r,hr)W{[l + z]E 7 r) 

= J dr EW([1 + z]E,r) 



J d 3 x J dM Y^ s h{x-Xi) 



M 

x S D (M - Mi)- u{nr - x\M) 

i LftiPc 

(B2) 

where 6(r, fir) is the overdensity. Therefore, we need 
to evaluate the ensemble average of the product of seed 
densities as follows: 



S%(xi - x i )S D (M 1 - Mi) $d(*2 - Xj)S D (M 2 - Mj)J2 S D (xa - x k )S D (M 3 - M k ) 

\ i 

dn dn dn 



j k 



dM x dM 2 dM 3 
+^ 3) {M 1 ,M 2l M 3 ,x 1 ,x 2l x 3 ) 



l + &>(M 1 ,M 2 ,x 1 ,x 2 )+e s 2> (M 2 ,M 3 ,x 2 ,x 3 )+e s 2, (M 1 ,M 3 ,x 1 ,x 3 ) 



dn 


dn 


dM x 


dM 2 


dn 


dn 


dM 2 


dM 3 


dn 


dn 



dMi dM 3 



1 + ^{M 1 ,M 2 ,x 1 ,x 2 )\ 6 3 D {x 3 - x 2 )S D {M 3 - M 2 ) 
l + (i 2) (M 2 ,M 3 ,x 2 ,x 3 )\ S 3 D (x! - x 3 )6 D {Mi - M 3 ) 
1 + ^ (M ll M 3 ,x ll x 3 )] 5 3 D (x 2 - Xl )6 D (M 2 - Mi) 



dn 

dM~! 



S 3 D (x 2 - Xl )5 D {M 2 - M 1 )8 3 D {x 3 - Xl )S D (M 3 - Mi). 



(B3) 
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The first term represents the three-halo contribution, i.e., 
three points considered Xx,x%, and 033 are in three dif- 
ferent halos; the second to fourth terms are the two-halo 
contribution, where two points are in one halo and the 
rest one point is in the other halo; the last term shows 
the one-halo contribution, where all three points are in 
the same halo. Now, our particular focus here is that 
one point, say x\, represents the blazar position, while 



the other two, x 2 and X3, the place of the dark matter 
annihilation. Since the latter effect is proportional to the 
density squared at one point, we have x 2 = x$. In this 
case, considering the fact that the halos are spatially ex- 
clusive, one can omit terms except for the second (2-halo 
term) and the last (1-halo term). 

We then evaluate the one-halo and two-halo terms. By 
substituting the above expression, we get 



E 2 (I B (E,n 1 )I D (E,h 2 )) lh 



= 4^ /^i / dr 2 EW([l + z 2 ]E,r 2 ) J dC $(£)^(A *i) (t^j) u 2 {x 2 - x x \M[C]), 



E 2 {I B {E 1 fi 1 )I D {E,h 2 )) 2h 



= E 2 (I B (E))(I D (E)} + -L J d 3 Xl J dr 2 EW([l + z 2 ]E,r 2 ) J d£ <1>{£)T E {£, Zl ) J dM, 
x J d 3 yi u 2 (x 2 -yi\M 1 )^ 2 \M[£},M u x l ,y 1 ), 



dn(Mx) ( M 



dM \n mP , 



where we used Eq. (6) of AK06 to reach {I B {E)){I D {E)) 
in the second expression. We again use small separation 
approximation, d 3 Xidr 2 /4ir — r\dr\dr 2 = r 2 drdrj, and 
get 

£ 2 (/ B (£,ni)/ D (£,n 2 )) 
= J drdr, r 2 EW{[\ + z]E, r) J dC $(£)F E (£, z) 



dM 



, dn(M') ( M' 



d 3 y 



dM \n m p, 

xu 2 (x 2 - y|M')d 2) (M[C],M', x llV ) 



(B4) 



(B5) 



(B6) 



M[£) 

^rnPc 



! (r;r + rdO\M[£]) 



To obtain the angular power spectrum, we Fourier- 
transform this expression. Using the similar relations to 
Eq. (|A12jl . we finally arrive at Eqs. (J29J) and ([30]) . 
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